Depinning and dynamics of AC driven vortex lattices in random media 
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We study the different dynamical regimes of a vortex lattice driven by AC forces in the presence of 
random pinning via numerical simulations. The behaviour of the different observables is charaterized 
as a function of the applied force amplitude for different frequencies. We discuss the inconveniences 
of using the mean velocity to identify the depinnig transition and we show that instead, the mean 
quadratic displacement of the lattice is the relevant magnitude to characterize different AC regimes. 
We discuss how the results depend on the initial configuration and we identify new hysteretic effects 
which are absent in the DC driven systems. 



I. INTRODUCTION 



A remarkable diversity of physical systems belong to 
the category of driven elastic manifolds moving over ran- 
dom landscapes. Extensively studied examples are abun- 
dant in the literature: moving vortex lattices in typ e II 
superconductors [ll-[T2j , sliding colloidal particles [13, 
or charge density wave systems 11511 . magnetic bubble ar- 
rays |16j] . driven Wigner crystals |17|.IT8| and stripe form- 
ing systems [l^. Among the topics that have recently 
received much attention is the nature of the depinning 
mechanisms that occur as the external DC driving force 
is increased beyond the critical force F^'-^ , and its rela- 
tion with proliferation of topological defects in the form 
of bounded or unbounded disclinations. The dependence 
of the dynamics with initial conditions has been exam- 
ined in detail and it was determined that for DC drives, 
memory of initial conditions is lost at depinning [2^ . A 
dynamical depinnig transition has been identified at a 
given force Fp'-^ . For larger forces topological defects 
heal and smectic linear flow is observed [1, . 

It should be noted however, that most of the research 
thrust has been specially focused on DC drives. In exper- 
iments, an AC field is often applied to order the vortex 
lattice (VL); the most accepted picture is that the AC 
field assists the system in an equilibration process, from 
a disordered metastable configuration to the equilibrium 
Bragg Glass phase, free of dislocations ^2^. However, a 
large amount of results, in both experiments 22, 23] and 
simulations [isl. f24l - l27j . show that, in some cases, an AC 
drive can disorder the VL, and that oscillatory dynamics 
plays an essential role. 

The oscillatory dynamics of the VL is in itself a broad 
field that is not completely understood. The main por- 
puse of this work is to provide a more comprehensive de- 
scription, using numerical simulations, of AC driven vor- 
tex lattices over a random distribution of pinning sites. 
We explore the eff'ect of the applied force amplitude for 



different frequencies, starting from different initial config- 
urations. The VL mean velocity, which can be dephased 
from the excitation, is no longer an adequate observable 
to detect depinning and we show that instead, the mean 
quadratic displacement of the lattice is the relevant mag- 
nitude to identify depinning and dynamical AC regimes. 
The depinning transition becomes a crossover as is de- 
scribed below. The pinned linear Campbell regime (28j 
is also identified and simulations are compared with an- 
alytical calculations. 

The paper is organized as follows: In Sec. HI] we de- 
scribe in detail the numerical simulations; starting with 
the model (jll A[) . following with the numerical procedures 
(|IIBp and giving the definition of the observables pi Cp . 
In Sec. |lll]we present and discuss the results and in Sec. 
IIVI we highlight the main conclusions. 



II. NUMERICAL SIMULATIONS 

A. The model 

In our simulations, we consider rigid vortices with 
coordinates in a two-dimensional rectangle of size x 
Ly that evolve according to the dynamics 



?7v,j 



(1) 



where Vj its velocity and rj the Bardeen-Stephen viscosity 
coefficient and is given by the sum of the vortex- vortex 
interaction, the pinning attraction and the Lorentz force 



(2) 



by, 



The vortex-vortex interaction per unit length is given 



E 



F™(r.-r,), 



(3) 
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where 
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F™(r. - r,) = -^fM^-I^j^)^.r (4) 

Here, (jjQ is the quantum of magnetic flux, A is the London 
penetration length and Ki the special Bessel function. 
The parameter f^^y is dimensionless and can be related 
to the stiffness of the vortex lattice (see Ref. [23 ). The 
Np pinning centers are supposed to be located at ran- 
dom positions R^ , and their interaction with vortices is 
modeled by 

F''f(r,) = g = F''f(r,-R,), (5) 

F'"'(r, - R,) = -F/e-^^^)'(r, - R,), (6) 

here, Fj (chosen from a Gaussian distribution with mean 
value PP = 0.2 and a standard deviation of O.IF^') and 
Tp tune the strength and range of the interaction. 

The Lorentz force per unit length is given by Fl = 
4>oJext X z where Jext is the external driving current den- 
sity and z is the versor perpendicular to the plane. 

Following Ref. [13 we measure lengths in units of A, 

forces (per unit length) in units of fo = gj^°^;j , time in 
units of to = rjX/ fo and frequency in units of ujq ~ 1/io- 
In our simulation, we will consider — 1600, — 40A, 
Ly = Np = 25Ny and rp = 0.2A. 

Concerning the numerical details, the equations of mo- 
tion are integrated using a standard Euler algorithm with 
step h = 0.04^0, and a hard cut-off A = 4A in calculat- 
ing the vortex- vortex force. When calculating physical 
observables we average over four realizations of disorder. 



C. Observables 

It is by now well known that the mobility of the VL 
is affected by the topology of the VL configuration . A 
common observable used to characterize the VL configu- 
ration is the density of lattice defects, n^, (i.e. vortices 
with 5 or 7 neighbors in the Delaunay triangulation) . 

The mobility itself is characterized by the vortex ve- 
locity. The instantaneous mean vortex velocity in the 
direction of the applied force at time t is 

i 

In the DC case, the v vs curves can be directly 
related with experimental current- voltage characteristics 
{V ~ I curves). 

In the AC case the relationship between velocity and 
force and the experimental V — I curve is more subtle, 
as the phase between the applied force and the velocity 
plays an essential role. 

We then define the half cycle mean velocity in the n*'* 
cycle as 

2 rnT+T/2 

Vn = 7f <t)dt (9) 

^ JnT 

Notice that the half cycle velocity takes into acount 
this phase factor. We estimate the mean phase factor ipn 
by applying a sinusoidal F^ (t) and calculating the phase 
between F^ and the first harmonic of the mean velocity 
v{t) in the n*^ cycle. 

Moreover, a quantity of interest is the average 
quadratic displacement defined as 



B. Procedures 

The AC drive is simulated with an external square AC 
force of the form 

pAC _ \ h if nT <t <nT +^ , . 

^ 1 -/l nT+^ <t<in + l)T ^'^ 

We study vortex lattice evolution as we vary the force 
in a slow ramping of Jl starting from different initial con- 
ditions. In each case, we leave the system evolving from 
the initial configuration with zero external force for 5 cy- 
cles. This free evolution creates or annihilates a low den- 
sity of defects and becomes our metastable initial state. 
A numerical realization will be obtained by ramping from 
/l to /l + A/i , applying Na cycles this new force, and 
then allowing the system to relax for Nu, cycles at zero 
force, before reassuming the ramping. Measuring is per- 
formed during the last cycle, just before switching off the 
applied force (i.e. during cycle Na). 



< SX^ > {t2,h) = ^ Y.^x,{t2) - x,{h)f. (10) 

In the AC case, we can take t2 — ti + pT, the mean 
quadratic displacement in the force direction after per- 
forming p oscillations. This quantity, as we will show 
below, is very useful for describing AC drives. 

III. RESULTS AND DISCUSSION 

Throughout the work we have explored the behaviour 
of the system as a function of the AC amplitude. Recent 
numerical studies on a similar system (colloids under the 
influence of an external AC drive and quenched disorder) 
[isj have found a crossover from a low drive-highly dis- 
orderd phase to a high drive-low disordered phase. In all 
cases considered in Ref. the initial configuration was 
highly disordered. Our results confirm the existence of 
the high drive-low disorderd phase but we will show that 
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FIG. 1: Number of topological defects (a) and half cycle mean 
velocity Vn (b) (see text) as a function of the amplitude of the ap- 
plied Lorentz force for different frequencies ui. All the curves are 
calculated slowly ramping f^, from an ordered VL configuration. 
Inset; zoom of HuC/l) curves at low amplitudes: > at finite 
frequency. 

when more general initial conditions are considered then 
richer, low drive phases emerge. 

Results starting from an ordered configuration are 
summarized in Figure [T] and Figure [51 In all the cases 
Na = Q cycles have been applied with each force. The 
observables w„ (Eq. [S]) and Ud have been calculated in 
the n"* cycle before switching off Jl- In Figure[T]we show 
the behaviour of the defect density Ud (Fig. [ija)) and the 
characteristic w„ — Jl curve (Fig. (Hb)) for a particular 
choice of vortex-vortex interaction ( fyy — O.8/0), at vari- 
ous frequencies uj. In the DC limit (shown in red full cir- 
cles) we clearly identify three different regimes, a pinned 
lattice for /l < F^'~^ , a disordered flow (plastic) region 
for F^^ < Jl < Fp'^ and a flowing (smectic) linear 
regime for Jl > F^'^' consistent with previous work . 
While the existence of F^'^' can be inferred form these 
AC curves, as the value for which all curves that present 
plastic motion merge, the distinction between the pinned 
and plastic regimes is less apparent. This is so, because 
even a weak driving force produces oscillations of vortices 
around their pinning centers with finite mean velocity. 
This is a strong argument to abandon the mean velocity 
as the adequate observable to indicate depinning. The 
inset of Fig. \Vp shows an enlarged area for small driving 
forces: at very small AC amplitudes a linear Campbell 



regime holds but the response becomes non-linear at am- 
plitudes well bellow F^^ . The definition of a critical 
force in the AC response is therefore not so obvious, and 
we will discuss this point in detail later. 

In the AC curves, we distinguish two frequency 
regimes: at low frequencies (w < 0.02wo in our simu- 
lation), restoring and pinning forces prevail over losses, 
whereas at high frequencies viscous drag dominates. The 
low frequency regime is characterized by a highly nonlin- 
ear response at intermediate AC amplitudes, associated 
with a pronounced peak in n^. This feature smears at 
higher frequencies, the density of topological defects de- 
creases, and the non-linearity is less apparent, going to a 
linear "ohmic" regime (blue down triangles) at very high 
frequency (w > 0.2ujq) where dissipative forces govern. 
Memory and history effects occur when pinning and elas- 
tic forces compete and both prevail over viscous forces. 
Therefore, in the following analysis we will focus in the 
low frequency regime. 




fL 



FIG. 2: Vn (a-left), phase factor (p (a-right), n^j (b-left) and mean 
quadratic displacement < SX^ > (b-right) as a function of fi^ at 
u] = O.OlSSoiQ. Vertical dashed lines indicate the different AC re- 
gions discussed in the text. The depinning of the VL occurs in 
a crossover region between Fd and Fc2, where < 5X'^ > grows. 
Inset: Comparison between riii before (black full circles) and after 
(grey full circles) N-w cycles without applied force (see text). Be- 
yond the depinning region the moving lattice configurations are 
unstable and relax. 

In Figure [5] we plot together various observables as 
a function of the AC amplitude for a fixed frequency 
uj — O.OISSwq, in a procedure analogous to that explained 
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in the description of Figure [T] Figure [2ja) displays the 
modulus of the half cycle mean velocity z;„ (left axis) and 
the estimated phase factor (right axis). In Figure 
[^fb) the density of defects rid (left axis) and the mean 
quadratic displacement < SX^ > after 5 cycles (right 
axis) are shown. The dashed horizontal line in Figure 
[5Jb) indicates de density of defects n^r that are sponta- 
neously created (without any applied force), and corre- 
spond to the initial more ordered configuration. A still 
more ordered configuration is unstable. 

We identify various regimes. At very low AC am- 
plitudes (/i < F;) a linear Vnifh) holds. The veloc- 
ity is mainly out of phase. The VL configuration re- 
mains unchanged, because vortices are trapped and can 
only perform small harmonic oscillations around their ini- 
tial positions. This motion is reversible, and therefore 
< 5X^ >— 0. In the other limit, at very high AC am- 
plitudes (/i > Fp), there is a dynamic reordering and 
the pinning potential is completely smeared. There is a 
linear in phase Ohmic response. 

In the intermediate range Fi < fi^ < Fp, the response 
is non-linear. A very rich behaviour with different non- 
linear regimes may be observed. In the first non-linear 
region, at small amplitudes { Fi < Jl < Fd), the veloc- 
ity is small, and there is not an appreciable displacement 
in the direction of the force (i.e. < 5X^ >~ 0); plastic 
random displacement produces a huge number of dislo- 
cations. However, most of the vortices remain trapped 
around the pinning sites. Depinning occurs in a crossover 
region, between Fd and Fc2- We identify the beginnig 
of the depinning region at the force Fd , where vortices 
move in average distances larger than the pinning radius 
Tp. In this region, the irreversible displacement grows 
with Jl and is reflected in the growth of < 5X^ >, that 
reaches its maximun at the second force Fc2 ; as this hap- 
pens, the density of defects attains its maximum and the 
slope oivn{fL) grows. Above Fc2, in the upper non- linear 
region {Fc2 < /l < -fp), the mean phase factor drasti- 
cally decreases (i.e. dissipation becomes relevant), while 
the motion becomes more reversible and the VL more or- 
dered. This is accompanied by a smooth decrease in the 
VnifL) slope that approaches the final linear relationship. 
Another interesting point to remark is that VL configura- 
tions that have similar density of defects at both sides of 
the depinning transition are qualitatively different. This 
difference can be detected, observing the relaxation of 
the VL configurations after removing the applied force. 
In the inset of Figure [2l^b), we compare the density of de- 
fects Ud before switching off the applied force (black full 
circles) with the same observable, cycles after remov- 
ing the applied force (grey full circles). While for forces 
lower than Fc2 , the vortices move around pinning sites in 
robust metastable configurations, beyond the depinning 
transition the moving lattice configurations are unstable, 
and relax towards a state with more defects. 

We now turn to analyze the various AC regimes start- 
ing from different initial conditions. In Figure[3]the den- 
sity of defects Ud (Fig. WL^)) and the mean quadratic 



^C1 ^C2 



(a) 






- 








fvv = 0.8fo 




^ □ 

tj 


n 


CO = 0.0188coo 








1 




1 ... 1 ... 1 . . 




0.016 0.024 
f. 



FIG. 3: (a) rid as a function of starting from different ini- 
tial configurations. Arrows indicate tfie direction of change of /i. 
Above Fc2 all the ridifL] curves taken with increasing merge. 
An hysteresis is observed (full and empty magenta squares), (b) 
< > (/l) curves corresponding to the same processes de- 

scribed in (a). This observable is independent of the initial condi- 
tion for the whole driving force range. 



displacement < SX"^ > (Fig. ^h)) as a function of the 
AC force are shown for different initial conditions and 
AC histories. The results described in the previous para- 
graph (starting from an ordered VL and increasing the 
AC force) are plotted in red crosses. Analogous results 
obtained with a similar AC protocol but starting from 
the more disordered metastable configuration are plot- 
ted in blue dots, and results with an initial intermediate 
Ud are plotted in green triangles. As in the DC case stud- 
ied in Ref. [l^, the different Ud curves merge for Jl larger 
than a threshold value. For the DC case, this threshold 
value was identified with I^'^ , the critical current at de- 
pinning. As we mentioned before, the identification of 
the critical current I^'^' is less obvious in the AC case, 
and we notice that the nd{fL) curves merge around Fc2- 
In magenta squares, we show the behaviour of Ud for 
decreasing (open squares) and a subsequent increasing 
(full squares) external driving. We observe the presence 
of hysteresis in the region between Fd and Fp . Below 
Fd pinning forces dominate, and above Fp they are com- 
pletely smeared. In the intermediate region, the strong 
competition gives rise to hysteresis. This hysteretic ef- 
fect beyond the critical force has not been observed in 
the DC casefl^ and it seems to be a new geninue effect 
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of the AC case. 

As another striking result, the behaviour of < > 
vs. Fl is independent of the initial condition for the 
whole driving force range. The reorganization of vortex 
defects by AC drives in initially different configurations, 
and for different AC protocols, surprisingly does not in- 
volve different quadratic displacements. This observable 
is univocally related with the applied force, and seems to 
be the best one to characterize the AC regimes. 
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FIG. 4: Normalized half cycle mean velocity Vn/ fh as a function 
of frequency ui in the linear regime, at low Fj^. Fit with EQ ll4l (con- 
tinuos line), considering a Gaussian distribution for the Labusch 
constant, gives ^ ~ a;c = 0.11 ± 0.05. 

Finally, we discuss in more detail the linear regime ob- 
tained at very small AC driving forces. In the harmonic 
approximation we can write the mean restoring force (av- 
eraging F™ -I- F"^ in Eq. [2]) as Frest= —ax where x 
denotes the average displacement of vortices from their 
equilibrium position and a is the Labusch parameter [28| . 
Within this aproximation, Eq. [1] can be written as: 

r]x + ax = F^^. (11) 

If we consider a semicycle for which F^'~' is positive, 
the solution of equation pT|) is 

x{t) + (xoiU,) - :^)e-t(*-*"), (12) 

a a 

where nT < t < nt + 

This equation provides a recurrence relation for xo{tn) 
that can be solved explcitly. In the limit n — > oo (the 
stationary regime) this gives 

a:o(tn) = -— tanh(^) (13) 

and the average velocity over half-cycle as 

2a; fr , , an . 
^= ^iitanh( — ). (14) 
an Zrjuj 

On the other hand, numercially we have performed dif- 
ferent calculatons of the half cycle mean velocity in the 



n*'' cycle f„ (according to Eq. |9]) as a function of u; for 
different values of /l, so we can average over and 
study the behaviour of V^J /l . In addition, as we did 
before, we have performed averages over different real- 
izations of pinning centers. This fact implies that indeed 
a is a random variable that as a first aproximation can 
be assumed to be Gaussian distributed and characterized 
by its mean value a and standard deviation a a ■ 

In Figure[4]we show the results of the simulation, where 
we have used 7 values of /l from 0.00025/o to 0.004/o 
and 4 realizations over disorder. Our fit gives a = 0.11, 
= 0.05. 

A critical frequency may be estimated as Wc ~ ^ = 
0.11 ± 0.05. We are able now to justify the two frequency 
regimes introduced at the beginning of this section. 

IV. CONCLUSIONS 

We have presented a comprehensive description of the 
dynamics of AC driven vortex lattices over a random dis- 
tribution of pinning sites, focusing the study in the low 
frequency regime, uj << ujc, where pinning and elastic 
forces prevail over viscous forces. This critical frequency 
ujc has been estimated by fitting results of numerical sim- 
ulations at low AC amplitudes with that predicted by an 
analytical model. 

The reorganization of vortex defects by AC drives from 
initially different configurations and for different AC pro- 
tocols produce history dependent vortex lattice configu- 
rations with different density of defects and mobility. 

As an important result, we notice that the behaviour 
of the mean quadratic displacements < SX"^ > (/l) is 
independent of the initial condition for the whole driving 
force range. This observable is univocally related with 
the applied force, and seems to be the best one to char- 
acterize the AC regimes. This should be contrasted with 
the mobility that can be large even in pinned AC driven 
lattices. 

As expected, a linear Campbell regime holds at very 
low AC amplitudes { Jl < Fi ) whereas the response is 
Ohmic at very high AC amplitudes (/^ > Fp ) where 
there is a dynamic reordering and the pinning potential is 
completely smeared. In both linear responses the motion 
is reversible and < 5X^ >= 0. 

In all the intermediate range Fi < Jl < Fp, the re- 
sponse is non-linear and a very rich behaviour with dif- 
ferent non-linear regimes may be observed. 

Depinning occurs in a crossover region {Fd < /l < 
Fc2) where < SX^ > grows with and vortices move in 
average distances much larger than the pinning radius. 

In the first non-linear region, below the depinning tran- 
sition { Fi < Jl < Fci) there is not an appreciable dis- 
placement in the direction of the force (i.e. < SX'^ 0); 
plastic random displacement produces a huge number 
of dislocations but most of the vortices remain trapped 
around the pinning sites. The density of defects Ud is 
strongly dependent on the initial configuration. 
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Beyond the depinning region the memory of the initial 
configuration is lost. Dynamics becomes more and more 
reversible as < 5X'^ > decreases. Dissipation becomes 
relevant and the VL reorders. The moving vortex lattice 
configurations become unstable and relax after switching 
off the applied force;. However, there is a large region ( 
Fc2 < Jl < Fp) where plastic motions are still present, 
as can be inferred from < 5X'^ >> and the non- linear 
response. Below the plastic threshold Fp, an hystere- 
sis in the applied force, absent in DC driven lattices, is 



observed, showing once again the rich particularities of 
oscillatory dynamics. 
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